Paradoxical diffusion: Discriminating between normal and anomalous random walks 



Bartlomiej DybiecQ and Ewa Gudowska-Nowal|3 

Marian Smoluchowski Institute of Physics, and Mark Kac Center for Complex Systems Research, 
Jagellonian University, ul. Reymonta 4, 30-059 Krakow, Poland 
(Dated: October 16, 2009) 

Commonly, normal diffusive behavior is characterized by a linear dependence of the second cen- 
tral moment on time, (x 2 (t)) oc t, while anomalous behavior is expected to show a different time 
dependence, (x 2 (t)} oc t s with S < 1 for subdiffusive and S > 1 for superdiffusive motions. Here we 
demonstrate that this kind of qualification, if applied straightforwardly, may be misleading: There 
are anomalous transport motions revealing perfectly "normal" diffusive character ({x^tt)) oc t), yet 
being non-Markov and non-Gaussian in nature. We use recently developed framework [j], Phys. Rev. 
E 75, 056702 (2007)] of Monte Carlo simulations which incorporates anomalous diffusion statistics 
in time and space and creates trajectories of such an extended random walk. For special choice 
of stability indices describing statistics of waiting times and jump lengths, the ensemble analysis 
of paradoxical diffusion is shown to hide temporal memory effects which can be properly detected 
only by examination of formal criteria of Markovianity (fulfillment of the Chapman-Kolmogorov 
equation) . 

PACS numbers: 05.40.Fb, 05.10.Gg, 02.50.-r, 02.50.Ey, 



I. INTRODUCTION 



Usually various types of diffusion processes are classi- 
fied by analysis of the spread of the distance traveled by 
a random walker. If the mean square displacement grows 
like ([x — x(0)] 2 ) oc t s with 5 < 1 the motion is called sub- 
diffusive, in contrast to normal (5 = 1) or superdiffusive 
(6 > 1) situations. For a free Brownian particle moving 
in one dimension, a stochastic random force entering its 
equation of motion is assumed to be composed of a large 
number of independent identical pulses. If they posses 
a finite variance, then by virtue of the standard Cen- 
tral Limit Theorem (CLT) the distribution of their sum 
follows the Gaussian statistics. However, as it has been 
proved by Levy and Khintchine, the CLT can be general- 
ized for independent, identically distributed (i.i.d) vari- 
ables characterized by non-finite variance or even non- 
finite mean value. With a Levy forcing characterized by 
a stability index a < 2 independent increments of the 



Darticle position sum up yielding ([x — x(0)Y 



x 



t 2/a 



[2|, see below. Such enhanced, fast superdiffusive motion 
is observed in various real situations when a testpar- 
ticle is able to perform unusually large jumps [1, U @. 
Levy flights have been documented to describe motion 
of fluorescent probes in living polymers, tracer particles 
in rotating flows and cooled atoms in laser fields. They 
serve also as a paradigm of efficient searching strategies 
@j 0i Si m social and environmental problems with some 
level of controversy [{|. 

In contrast, transport in porous, fractal-like media or 
relaxation kinetics in inhomogeneous materials are usu- 
ally ultraslow, i.e. subdiffusive [H, E3, EH- The most 
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intriguing situations take place however, when both ef- 
fects - occurrence of long jumps and long waiting times 
for the next step - are incorporated in the same scenario 
[l2| . The approach to this kind of anomalous motion 
is provided by continuous time random walks (CTRW) 
which assume that the steps of the walker occur at ran- 
dom times generated by a renewal process. In particular, 
a mathematical idealization of a free Brownian motion 
(Wiener process W(t)) can be then derived as a limit (in 
distribution) of i.i.d random (Gaussian) jumps taken at 
infinitcsimally short time intervals of non-random length. 
Other generalizations are also possible, e.g. W(t) can 
be defined limit of random Gaussian jumps per- 

formed at random Poissonian times. The characteris- 
tic feature of the Gaussian Wiener process is the conti- 
nuity of its sample paths. In other words, realizations 
(trajectories) of the Wiener process are continuous (al- 
though nowhere differentiable) [l3[. The process is also 
self-similar (scale invariant) which means that by rescal- 
ing t' = Xt and W'(t) = A _1/2 W(Ai) another Wiener 
process with the same properties is obtained. Among 
scale invariant stable processes, the Wiener process is the 
only one which possesses finite variance [lj, E3, EE EH • 
Moreover, since the correlation function of increments 
AW(s) = W(t + s) — W(t) depends only on time differ- 
ence s and increments of non-overlapping times are sta- 
tistically independent, the formal differentiation of W(t) 
yields a white, memoryless Gaussian process [l7| : 



(i) 



commonly used as a source of idealized environmental 
noises within the Langevin description 



dX{t) = f(X)dt + dW(t). 



(2) 
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Here f(X) stands for the drift term which in the case of 
a one-dimensional ovcrdamped motion is directly related 
to the potential V(X), i.e. f(X) = -dV(X)/dX. 
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In more general terms the CTRW concept may asymp- 
totically lead to non-Markov, space-time fractional noise 
and in effect, to space-time fractional diffusion. For 
example, let us define AW(t) = AX(t) = J^fJi X n 
where the number of summands N(t) is statistically 
independent from Xi and governed by a renewal pro- 
cess EfJi T i < t < Efi? 41 ^ with t > 0. Let 
us assume further that Tj, Xi belong to the domain 
of attraction of stable distributions, T, ~ S v ,i and 
Xi ~ <Sa,/3, whose corresponding characteristic functions 
<j){k) = (exp(ikS ai p)) = f x oo e lhx l a .p(x](j)dx, with the 
density l a ,p(%] cr), are given by 



where c\ and C2 are constants. The resulting (in gen- 
eral non-Markov) process becomes y/a self-similar Levy 
random walk 0, M, E, EI H HH, Ha, H3 , i.e. 



P (x,t) = r"/>(ir"/ ,i). 



(9) 



The asymptotic form given by Eq. (jSJ) can be easily 
derived [ll, M, HI for decoupled CTRW by applying 
Tauberian theorems to the Montroll- Weiss [26[ expres- 
sion 



p(q, u) 
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for the Laplace- Fourier transform of p(x,t). The latter 
satisfies the integral (master) equation 



for a = 1. Here the parameter a £ (0,2] denotes the 
stability index, yielding the asymptotic long tail power 
law for the x-distribution, which for a < 2 is of the 
|a;|~( 1+Q type. The parameter a (a £ (0, oo)) charac- 
terizes the scale whereas j3 ((3 £ [—1,1]) defines an asym- 
metry (skewness) of the distribution. 

Note, that for < v < 1, f3 = 1, the stable variable 
S^i is defined on positive semi-axis. Within the above 
formulation the counting process N(t) satisfies 

r Arm i f [{t,c) " x] 

lim Prob { , // < x \ = lim Prob < V T t > t 

\(t/cY J t-oc ^ 
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Here ij)(t) stands for a waiting time probability density 
function (PDF) whereas w(x) denotes a jump length 
PDF. With a suitable change of time and space vari- 
ables, in the limit of x — > oo, t — > oo, the Laplace-Fourier 
transform p{q 1 u) can be written as 



p{q,u) 



lim Prob < V T t > — 

n— >oc I ^ — ' x 



(5) 



dsexp[-s(u v + \q\ a )}. (12) 
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The inverse Laplace transform of p(g, u) can be expressed 
in a series form fijl : 



where [{t/c) v x\ denotes the integer part of the number 
{t/c) v x and L a ^(x) stands for the stable distribution of 
random variable S a< p, i.e. l a ,p{x) = dL a ^{x)/dx. More- 
over, since 



= g (-1)" 



ut du 



k=0 
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lim Prob 



L a ,p(x) (6) 



= E u (-\q\ a n, (13) 
where E v {z) is the Mittag-Leffler function of order v 



and 



k=0 



T(kv + 1) 



(14) 



P(x,t) = ^2p{x\n)p n {n(t)), 

n 

asymptotically one gets 



(7) 



Further application of the inverse Fourier transform in q 
yields [15| a final series representation of p(x, t): 



P (X, t) ~ (C 2 t)-^ a j l a> p ({c 2 t)-»/ a XT»/ a ) l uA (T)T"/ a dT, 

(8) 



1 ~ (-l) fc r(fcq + l) 
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where y = xjt v l a . The above series are divergent for 
a ^ v. However, for a = v, the summation has been 
shown [l5[ to produce the closed analytical formula 



p(x,t) 



sin(7r^/2) 



Av\t |y|" + |2/|-' y + 2cos(^/2) 



(16) 



where (as previously) y = x/t v ' a . 

The exact solution of the decoupled CTRW, as given 
by the infinite sum ([7|) of stable probability densities, has 
been studied by Barkai [2l[ based on a special choice of 
PDFs w(x) and ip(t). In particular, in [2l[ the extremely 
slow convergence of certain CTRW solutions to the (frac- 
tional) diffusion approximation has been discussed. The 
rigorous proof of equivalence between some classes of 
CTRWs and fractional diffusion has been given by Hilfer 
and Anton [27[ 

In this article we investigate CTRW scenarios which, 
in an asymptotic limit, yield paradoxical diffusion, i.e. 
the non-Markovian supcrdiffusive process taking place 
under sublinear operational time. The combination of 
long flights and long breaks between them is responsi- 
ble for the characteristic shape of the PDF and scaling 
properties of moments. In particular, for a = 2is, the 
paradoxical diffusion process exhibits the same scaling 
as an ordinary Brownian motion despite its PDF is sig- 
nificantly different from Gaussian. 

In a forthcoming Section, the relation between frac- 
tional calculus and CTRW approach is briefly reminded 
and the experimental results based on numerical PDF 
estimators are presented. Section III is devoted to the 
discussion of scaling properties of moments. In Section 
IV detection and analysis of memory effects in empirical 
series of the CTRW-type realizations are proposed and 
critically tested. 



II. RELATION BETWEEN CTRW AND 
FRACTIONAL CALCULUS 

The theory of stochastic integration of a corresponding 
Ito-Langevin equation with respect to a general CTRW 
"measure" dW(t) = dX has been developed in a series of 
papers 0, [H, HH Hil ■ Here, we study statistical prop- 
erties of such a motion constrained to the initial posi- 
tion X(0) = 0. To achieve the goals, we adhere to the 
scheme of stochastic subordination [l], [2iJ [3(j, i.e. we 
obtain the process of primary interest X(t) as a function 
X(t) = X(St) by randomizing the time clock of the pro- 
cess X(s) using a different clock St. In this approach St 
stands for ^-stable subordinator, St = inf {s : U(s) > t}, 
where U(s) denotes a strictly increasing z^-stable pro- 
cess whose distribution L Vi \ yields a Laplace transform 
^ e -kU(s)^ _ g-sfc 1 ^ rpjy, p arcn t process X(s) is composed 
of increments of symmetric a-stablc motion described in 
an operational time s 



and in every jump moment the relation U(St) = t is ful- 
filled. The (inverse-time) subordinator St is (in general) 
non-Markovian hence, as it will be shown, the diffusion 
process X(St) possesses also some degree of memory. The 
above setup has been recently proved [l], H^, [3(| to give 
a proper stochastic realization of the random process de- 
scribed otherwise by a fractional diffusion equation: 



dp(x, t) 



d d a 



p(x,t), (18) 



with the initial condition p(x,0) = S(x). In the above 
equation o-Df 1- ^ denotes the Riemannn-Liouville frac- 
tional derivative qD\~ v = ■§tqD^' / defined by the relation 



D\- V f{x,t) = 



1 d 



Y{v) dt J (t - t> ) 



(19) 



and 



— stands for the Riesz fractional derivative 

with the Fourier transform ^F[ d ^ x ^ f( x )] = ~ \k\ a f(x). 
Eq. (Tl8|) has been otherwise derived from a generalized 
Master equation [311] . The formal solution to Eq. 
can be written 13111 as: 



p(x, t) = E v 



d , d a 



t»)p(x,0). (20) 



For processes with survival function ^(t) = 1— f Q ip(r)dT 
(cf. Eq. dTTJ) ) given by the Mittag-LefHer function 
Eq . (fljll) , this solution takes an explicit form [l5|, [24|, [2f| 



(21) 



n=0 



dX(s) = -V'(X(s))ds + dL afi (s), 



(17) 



where Eu(z) = 4^E v (z) and w n (x) cx l a o(x), see [24|, 

In this paper, instead of investigating properties of an 
analytical solution to E q. (1181) . we switch to a Monte 
Carlo method [l|, [29|, 134l.l35l] which allows generating 
trajectories of the subordinated process X(t) with the 
parent process X(s) in the potential free case, i.e. for 
V(x) = 0. The assumed algorithm provides means to ex- 
amine the competition between subdiffusion (controlled 
by a ^-parameter) and Levy flights characterized by a 
stability index a. From the ensemble of simulated trajec- 
tories the estimator of the density p(x, t) is reconstructed 
and statistical qualifiers (like quantiles) are derived and 
analyzed. 

As mentioned, the studied process is vja self-similar 
(cf. Eq. ^j). We further focus on examination of a spe- 
cial case for which vja = 1/2. As an exemplary values 
of model parameters we choose v = 1, a = 2 (Markovian 
Brownian diffusion) and v = 0.8, a = 1.6 (subordination 
of non-Markovian sub-diffusion with Levy flights). Ad- 
ditionally we use v — 1, a = 1.6 and v — 0.8, a = 2 
as Markovian and non-Markovian counterparts of main 
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cases analyzed. Fig. Q] compares trajectories for all ex- 
emplary values of v and a. Straight horizontal lines (for 
v = 0.8) correspond to particle trapping while straight 
vertical lines (for a = 1.6) correspond to Levy flights. 
The straight lines manifest anomalous character of diffu- 
sive process. 

To further verify correctness of the implemented ver- 
sion of the subordination algorithm we have per- 
formed extensive numerical tests. In particular, some 
of the estimated probability densities have been com- 
pared with their analytical representations and the per- 
fect match between numerical data and analytical results 
have been found. Fig. [5] displays numerical estimators of 
PDFs and analytical results for v = 1 with a = 2 (Gaus- 
sian case, left top panel), v = 1 with a = 1 (Cauchy 
case, right top panel), v = 1/2 with a = 1 (left bottom 
panel) and v = 2/3 with a = 2 (right bottom panel). For 
those last two cases, the expressions for p(x, t) has been 
derived [l5[ , starting from the series representation given 
by Eq. (fT5|) . For v = 1/2, a = 1 the appropriate formula 
reads 



KM) = -2^ eXP (l' 
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FIG. 2: (Color online) PDFs for v = 1 with a = 2 (left top 
panel), i/ = 1 with a = 1 (right top panel), v = 1/2 with 
a = 1 (left bottom panel) and ^ = 2/3 with a = 2 (right 
bottom panel). Eq. 1)171) was numerically approximated by 
subordination techniques with the time step of integration 
At = 10 -3 and averaged over iV = 10 6 realizations. Solid 
lines present theoretical densities: Gaussian (left top panel), 
Cauchy (right top panel) and the p(x,t) given by Eqs. (|22[) 
(left bottom panel) and Q23p (right bottom panel). Note the 
semi-log scale in the bottom panels. 



while for ^=-2/3, a ~2 the probability density is 

3 2/3 



Ai 



(3t)V3 



(23) 



Ei(x) and Ai(x) are the integral exponential function and 
the Airy function respectively. We have also compared 
results of simulations and Eq. (fTtjj) for other sets of pa- 
rameters v, a. Also there, the excellent agreement has 
been detected (results not shown). 
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FIG. 3: (Color online) PDFs (top panel) and 1 - CDF(a;, t) 
(bottom panel) at t = 2 (left panel), t — 15 (right panel). 
Simulation parameters as in Fig. [1] Eq. (|17p was numerically 
approximated by subordination techniques with the time step 
of integration At = 1CP 3 and averaged over iV = 10 6 realiza- 
tions. Solid lines present theoretical asymptotic x~ is scaling 
representative for a = 1.6 and v = 1, i.e. for Markovian Levy 
flight. 



FIG. 1: Sample trajectories for v = 1, a = 2 (left top panel), 
v = 1, a = 1.6 (left bottom panel), v — 0.8, a = 2 (right top 
panel) and v = 0.8, a = 1.6 (right bottom panel). Eq. (|17[) 
was numerically approximated by subordination techniques 
with the time step of integration At — 10~ 2 and averaged 
over = 10 6 realizations. 



Figure[3]and[4]display time-dependent probability den- 
sities p(x,t) and corresponding cumulative distribution 
functions (CDF(x,t) = /* 00 p(ar' ,t)dx') for "short" and 
for, approximately, an order of magnitude "longer" times. 
The persistent cusp (3(| located at x — is a finger-print 
of the initial condition p(x,0) = 5(x) and is typically 
recorded for subdiffusion induced by the subordinator St 
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FIG. 4: (Color online) PDFs (top panel) and 1 - CDF(a:,t) 
(bottom panel) at t = 2 (left panel), t = 15 (right panel). 
Simulation parameters as in Fig. [1] Eq. (|17p was numerically 
approximated by subordination techniques with the time step 
of integration At = 1CP 2 and averaged over N = 10 6 realiza- 
tions. Solid lines present theoretical asymptotic x~ 1,2 scaling 
representative for a — 1.2 and v = 1, i.e. for Markovian Levy 
flight. 



with v < 1. For Markov Levy- Wiener process [53, HH for 
which the characteristic exponent v = 1, the cusp disap- 
pears and PDFs of the process X(St) become smooth at 
x — 0. In particular, for the Markovian Gaussian case 
[y = 1, a = 2) corresponding to a standard Wiener diffu- 
sion, PDF perfectly coincides with the analytical normal 
density N(0,y/i). 

The presence of Levy nights is also well visible in 
the power-law asymptotic of CDF, see bottom panels of 
Figs. [3] and 0] Indeed, for a < 2 independently of the 
actual value of the subdiffusion parameter v and at arbi- 
trary time, p(x,t) oc for x — > oo. Furthermore, 
all PDFs are symmetric with median and modal values 
located at the origin. 



III. SCALING PROPERTIES OF MOMENTS 

The vja self-similar character of the process (cf. 
Eq. ©) is an outcome of allowed long flights and long 
breaks between successive steps. In consequence, the 
whole distribution scales as a function of x/t v ' a with 
the width of the distribution growing superdiffusively for 
a < 2v and subdiffusivcly for a > 2v. This W Q scal- 
ing is also clearly observable in the behavior of the stan- 
dard deviation and quantiles q p (t), defined via the rela- 
tion Prob{A(i) $C q P (t)} = p, see Figs. [51151 For random 
walks subject to supcrdiffusive, long-ranging trajectories 
(a = 1.6), the asymptotic scaling is observed for suffi- 
ciently long times, cf. Fig. [5] On the other hand, normal 
(Gaussian) distribution of jumps superimposed on sub- 
diffusive motion of trapped particles (y = 0.8) clearly 
shows rapid convergence to the vja law. Notably, both 
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FIG. 5: (Color online) Time dependence of vax(x)/t. Straight 
lines present t 2 "/ a_1 theoretical scaling (see Eq. (|30|l and ex- 
planation in the text). Simulation parameters as in Fig. [T] 



sets v = 1, a = 2 and v = 0.8, a = 1.6 lead to the same 
scaling i 1 / 2 , although in the case v = 0.8, a = 1.6, the 
process X(t) = X(S t ) is non-Markov, in contrast to a 
standard Gaussian diffusion obtained for v = l,a = 2. 
Thus, the competition between subdiffusion and Levy 
flights questions standard ways of discrimination between 
normal (Markov, (\x — x(0)] 2 ) oc t) and anomalous (gen- 
erally, non-Markov ([x — x(0)] 2 ) oc t s ) diffusion processes. 

Indeed, for v = 1, the process X(t) is not only 1/a 
self-similar but it is also memory less (i.e. Markovian). 
In such a case, the asymptotic PDF p(x, t) is ev-stable 
[3, [37], [H, H(| with the scale parameter a growing with 
time like t x l a . cf. Eq. |3]). This is no longer true for 
subordination with v < 1 when the underlying process 
becomes non-Markovian and the spread of the distribu- 
tion follows the ^/"-scaling (cf. Fig. [6J right panels). 

Some additional care should be taken when discussing 
the scaling character of moments of p(x,t) 0, H3, ItH ]. 
Clearly, Levy distribution (with a < 2) of jump lengths 
leads to infinite second moment (see Eqs. ([3]) and Q) 

/oo 
x l a ,p(x\ &)dx = 00, (24) 
-00 

irrespectively of time t. Moreover, the mean value (x) 
of stable variables is finite for a > 1 only ((x) =0 for 
symmetric case under investigation) . Those observations 
seem to contradict demonstration of the scaling visible in 
Fig. [S] where standard deviations derived from ensembles 
of trajectories are finite and grow in time according to a 
power law. A nice explanation of this behavior can be 
given following argumentation of Bouchaud and Georges 
[2|: Every finite but otherwise arbitrarily long trajectory 
of a Levy flight, i.e. the stochastic process underlying 
Eq. (|18|) with v = 1, is a sum of finite number of inde- 
pendent stable random variables. Among all summed N 
stable random variables there is the largest one, let say 
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l c (N). The asymptotic form of a stable densities 

/ Q , /3 (x;a)oc.T-( 1+Q ), 



(25) 



together with the estimate for l c (N) allow one to estimate 
how standard deviations grows with a number of jumps 
N. In fact, the largest value l c {N) can be estimated from 
the condition 



N j l a> p(x)dx 
h c {N) 



1. 



(26) 



which locates most of the "probability mass" in events 
not exceeding the step length l c (otherwise, the relation 
states that l c {N) occurred at most once in N trials, 0]). 
Alternatively, l c {N) can be estimated as a value which 
maximizes probability that the largest number chosen in 
N trials is L 




N-l 
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FIG. 6: (Color online) Quantiles: qo.9, go. 8, <?o.7, Qo.6 (from top 
to bottom) for v = l,a = 2 (left top panel), v = 1, a. = 1.6 
(left bottom panel), v = 0.8, a — 2 (right top panel) and 
" — 1 v = 0.8, a = 1.6 (right bottom panel). The straight line 
presents theoretical t"^ a scaling. Simulation parameters as in 

Fig.m 



(27) 

By use of Eqs. (|26|) and (f25|) . simple integration leads to 



(28) 



l c (N) cx N 1/a . 



Due to finite, but otherwise arbitrarily large number of 
trials N, the effective distributions becomes restricted to 
the finite domain which size is controlled by Eq. (j28|) . 
Using the estimated threshold, see Eq. (|28|) and asymp- 
totic form of stable densities, see Eq. (|23|) . it is possible 
to derive an estimate of (x 2 ) 



x 2 l a , p (x)dx « (N^ a ) = N 2 ' 



Finally, after N jumps 

( X 2 ) N =N(x 2 ) <xiV 2/c \ 



(29) 



(30) 



Consequently, for Levy flights standard deviation grows 
like a power law with the number of jumps N. In our 
CTRW scenario incorporating competition between long 
rests and long jumps, the number of jumps N = N(t) 
grows sublinearly in time, N cx t v , leading effectively 
to (a; 2 ) at cx t 2v l a with < v < 1 and < a < 2. 
Since in any experimental realization tails of the Levy 
distributions are almost inaccessible and therefore effec- 
tively truncated, analyzed sample trajectories follow the 
pattern of the W Q scaling, which is well documented in 

Fig.m 



set v/a = 1/2, (e.g. v = 0.8, a = 1.6) when contrary to 
what might be expected, ([x— x(0)] 2 ) cx t, similarly to the 
standard, Markov, Gaussian case. This observation im- 
plies a different problem to be brought about: Given an 
experimental set of data, say time series representative 
for a given process, how can one interpret its statistical 
properties and conclude about anomalous (subdiffusive) 
character of underlying kinetics? The similar question 
has been carried out in a series of papers discussing use 
of transport coefficients in systems exhibiting weak er- 
godicity breaking (see [42j and references therein). 

To further elucidate the nature of simulated data sets 
for vjoi = 1/2, we have adhered to and tested formal 
criteria [It], EH, SI, EH defining the Markov process. 
The standard formalism of space- and time- continuous 
Markov processes requires fulfillment of the Chapman- 
Kolmogorov equation (ti > t 2 > t 3 ) 

P(xi,ti\x 3 ,t 3 ) = P(x 1 ,t 1 \x 2 , t 2 )P(x2,t 2 \x3, h)- 

(31) 

along with the constraint for conditional probabilities 
which for the "memoryless" process should not depend 
on its history. In particular, for a hierarchy of times 
h > <2 > t 3 , the following relation has to be satisfied 



IV. DISCRIMINATING MEMORY EFFECTS 



P(x 1 ,t 1 |x 2 ,t 2 ) = P(x 1 ,t 1 \x 2 ,t 2 ,x 3 ,t 3 ) 



(32) 



Clearly, by construction, for v < 1 the limiting 
"anomalous diffusion" process X(t) is non-Markov. This 
feature is however non-transparent when discussing sta- 
tistical properties of the process by analyzing ensemble- 
averaged mean-square displacement for the parameters 



Eqs. (f3Tj) and (f3"2"]) have been used to directly verify 
whether the process under consideration is of the Marko- 
vian or non-Markovian type. From Eq. (|31[) squared 
cumulative deviation Q 2 between LHS and RHS of the 
Chapman-Kolmogorov relation summed over initial (x 3 ) 
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and final (xi) states has been calculated [4J 



E 



P(xi,ti|x 3 ,t 3 ) 

- / ] P(xx, ii|x 2 , t 2 )P(x 2 , h\x 3 , t 3 ) 



(33) 



The same procedure can be applied to Eq. (|32|) leading 
to 



M 2 = Yj [p(xiM\x2M) 



-P(xi 1 t 1 \x 2 ,t 2 ,x 3 ,t 3 ) . (34) 

Figure [7| presents evaluation of Q 2 (top panel) and M 2 
(bottom panel) for t\ = 27 and t 3 = 6 as a function of 
the intermediate time t 2 = {7, 8, 9, ... , 25, 26}. It is seen 
that deviations from the Chapman-Kolmogorov identity 
are well registered for processes with long rests when sub- 
diffusion wins competition with Levy flights at the level 
of sample paths. The tests based on Q 2 (see Eq. ([33|) ) 
and M 2 (see Eq. (|34p) have comparative character. The 
deviations Q 2 and M 2 are about three order of magni- 
tudes higher for the parameter sets v = 0.8, a = 2.0 and 
v = 0.8, a = 1.6 than Q 2 and M 2 values for the Marko- 
vian counterparts with v = 1 and a = 2, a = 1.6, re- 
spectively. Performed analysis clearly demonstrates non- 
Markovian character of the limiting diffusion process for 
v < 1 and the findings indicate that scaling of PDF, 
p(x,t) = t~ 1 / 2 p(xt~ 1 / 2 , 1) or, in consequence, scaling 
of the variance ([x — x(0)] 2 ) oc t and interquantile dis- 
tances (see Fig. [5]) do not discriminate satisfactory be- 
tween "normal" and "anomalous" diffusive motions 
In fact, linear in time spread of the second moment does 
not necessarily characterize normal diffusion process. In- 
stead, it can be an outcome of a special interplay between 
subdiffusion and Levy flights combined in the subordina- 
tion X(t) — X(St). The competition between both pro- 
cesses is better displayed in analyzed sample trajectories 
X(t) where combination of long jumps and long trapping 
times can be detected, see Fig. Q] 



V. CONCLUSIONS 

Summarizing, by implementing Monte Carlo simula- 
tions which allow visualization of stochastic trajectories 
subjected to subdiffusion (via time-subordination) and 
supcrdiffusive Levy flights (via extremely long jumps in 
space), we have demonstrated that the standard mea- 
sure used to discriminate between anomalous and nor- 
mal behavior cannot be applied straightforwardly. The 
mean square displacement alone, as derived from the (fi- 
nite) set of time-series data does not provide full infor- 
mation about the underlying system dynamics. In order 
to get proper insight into the character of the motion, 
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FIG. 7: (Color online) Squared sum of deviations Q 2 , see 
Eq. J33j, (top panel) and M 2 , see Eq. ((34j) , (bottom panel) 
for ti = 27, tz — 6 as a function of the intermediate time 
ti. 2D histograms were created on the [—10, 10] 2 domain. 
3D histograms were created on the [—10, 10] 3 domain. Due 
to non-stationary character of the studied process the analy- 
sis is performed for the series of increments x(t + 1) — x(t). 
Simulation parameters as in Fig. [T] 



it is necessary to perform analysis of individual trajec- 
tories. Subordination which describes a transformation 
between a physical time and an operational time of the 
system [l], [46[ is responsible for unusual statistical prop- 
erties of waiting times between subsequent steps of the 
motion. In turn, Levy flights are registered in instan- 
taneous long jumps performed by a walker. Super- or 
sub- linear character of the motion in physical time is 
dictated by a coarse-graining procedure, in which frac- 
tional time derivative with the index v combines with 
a fractional spatial derivative with the index a. Such 
situations may occur in motion on random potential sur- 
faces where the presence of vacancies and other defects 
introduces both - spatial and temporal disorder (4?J ■ We 
believe that the issue of the interplay of super- and sub- 
diffusion with a crossover resulting in a pseudo-normal 
paradoxical diffusion may be of special interest in the 
context of e.g. facilitated target location of proteins on 
folding heteropolymers [48| or in analysis of single parti- 
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clc tracking experiments [42|, H^, \E$, [HTj] , where the hid- 
den subdiffusion process can be masked and appear as a 
normal diffusion. 
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